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. ^ ■ A method of a non-stationary description of tunneling of a particle through 

the one-dimensional and spherically symmetric rectangular barriers on the basis of 
analisis of multiple internal reflections of wave packets in relation on the barrier 
boundaries, named as the method of multiple internal reflections, is presented at 
►^ , the first time. For the one-dimensional problem the applicability of this method is 

ly^ \ proved, its specific features are analyzed. For the spherically symmetric problem 

cn ■ the amplitudes of the transmitted and reflected wave packets, times of tunneling 

^^ . and reflection in relation to the barrier are calculated using this method. The effect 

^^ ' of Hartman-Fletcher is analyzed. 
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C ■ effect of Hartman-Fletcher 

>: 

^j : 1 Introduction 

As the further development of the time analysis of tunneling processes considered in 
Pl El E] , here we present the non-stationary solution method for the problem of scattering 
of nonrelativistic particle on the spherically symmetric field which has rectangular barrier. 
Such approach uses the account of multiple internal reflections of fluxes in relation to the 
boundaries of the barrier (this approach further we name as the method of multiple internal 
reflections and we present such a method at the flrst time). 

An important speciflc feature of this method is a consideration of particle propagation 
on the basis of wave packets (WPs). Due to this one can fulflU the time analysis of particle 
propagation (or tunneling) and study in details these processes in interesting time moment 
or in relation to the concrete point of space. In result, the validation of method on the 
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basis of the account of the multiple internal reflections of WP is given correctly both for 
above-barrier and sub-barrier regions (unlike stationary consideration of such approach 
to the problem solution [HE]). 

At the beggining we consider the problem of tunneling of particle through an one- 
dimensional rectangular barrier. This problem is test one and allows to analyze the 
specific features of this method. The time analysis of the method of multiple internal 
reflections is fulfilled for the first time. 

The interesting perspective of this method is enough easy calculation (in comparison 
with the known stationary methods) of stationary wave functions (WFs) for the one- 
dimensional and spherically symmetric problems with barriers of more composite form 
than rectangular. It is shown weakly in the one-dimensional problem with one rectangular 
barrier and is shown more brightly in the problems with barriers of more composite 
forms. As an example we consider the one-dimensional problem with the double-humb 
rectangular barrier. 

Further the method is used for solving the problem of scattering of particle on the 
spherically symmetric field, which radial part has a rectangular barrier. For it the ampli- 
tudes of transmitted and reflected WPs, total times of tunneling and reflection in relation 
to the barrier are found. The effect of Hartman - Fletcher is considered. The expression 
for the S-matrix is presented in the form of the sum of two components corresponding to 
the amplitudes of transmitted and reflected WPs. Spherically symmetric problem with 
use of the method of multiple internal reflections is solved for the first time. 

2 The tunneling evolution of particle through the 
ID- rectangular barrier 

Let's consider the problem of particle moving along the x-direction and tunneling through 
a rectangular potential barrier (see Fig. P). We label the region I for x < 0, region II 
for < X < a and region III for x > a, accordingly. Solving the stationary Schrodinger 
equations, one can find the general solution for WF: 

for X < 0; 
ip{x) = \ ae~^^ + /3e^^, for < x < a; (1) 

for X > a; 




where k = ^y2mE, ^ = ^J2m{Vi — E), E and m are the energy and mass of a parti- 
cle, accordingly. The coefficients At, Ar, a and j3 can be calculated analytically, using 
requirements of the continuity of WF ^{x) and its derivative on each boundaries of the 
barrier. 

The tunneling evolution of particle can be described using non-stationary consideration 
of the propagating WP constructed on the basis of solution of the stationary Schrodinger 
equations: 

+ 00 

^(x, t)= j g{E - E) * v5(x, k)e-'^'/^dE, (2) 



where the weight amplitude g{E — E) satisfies to the requirement of the normalization 

/ \g{E — E)\^dE = 1, value E is the average energy. 



By inserting in the integral (2) instead of the total wave ^{x, k) the incident ipinc{x, k), 
transmitted iptr{x,k) or reflected part of WF ipref{x,k), defined by expression (1), we 
obtain the incident, transmitted or reflected WP, respectively. Considering only sub- 
barrier processes, we exclude the component of waves for above-barrier energies, including 
the additional transformation: 

g^E -E)^ G{E -E)= g{E - E)9{Vi - E). (3) 

The method of multiple internal reflections considers the process of WP propagation 
sequentially on steps of its propagation in relation to each boundary of the barrier. 

At the first step we consider WP in region I, which is incident upon the first (initial) 
boundary of the barrier. Let's assume, that this packet transforms into the WP, transmit- 
ted through this boundary and tunneling further in region II, and into the WP, reflected 
from boundary and propagating back in region I. Thus from the reasons of elementary 
causality we consider, that WP, tunneling in region II, is not reached the second (final) 
boundary of the barrier because of a terminating velocity of its propagation, and conse- 
quently at this step we consider only two regions I and II. Because of the same physical 
reasons to write the expression for this packet, we consider, that its amplitude should 
decrease in positive x-direction. For this we use only one item /5*exp(— ^x) in expression 
(1), throwing the second increasing item a * exp(^x) (in the opposite case we break the 
requirement of the finiteness of WF for an indefinitely wide barrier). In result, for the 
incident, transmitted and reflected WPs in relation to the first boundary one can write: 

+ 00 _ 

^i^^(x, t) = j G{E -E)* e^^^-^-s*ArfE, X < 0, 



iPl(x, t) = J G{E -E)* /5°e-«^-^-^*/^t/E < x < a, (4) 



+ 00 _ 

7/,! (x,t) = / G{E -E)* A%e-'''''-'^^/^dE, x <0. 



We assume, that during a particular interval of time considered at the first step, the 
transmitted WP tunnels in the sub-barrier region in the positive x-direction. In result, 
at this step we receive the flux for this packet, not equal to zero. Considering only the 
stationary expressions for WF for this step, we receive the zero fluxes in the barrier region 

Constructing from (4) total WF in regions I and II and using the condition of continuity 
of this WF and its derivative, we find unknown coefficients /3° and A^^. 

At the second step we consider WP, tunneling in region II and incident upon the 
second boundary of the barrier at point x = a. It transforms into the WP, transmitted 
through this boundary and propagated in region III, and into the WP, reflected from 
boundary and tunneled back in region II. Let's define them so: 

'^inc{x,t) = ipl{x,t), 0<x<a, 

+ 00 _ 

^Pl{x,t) = j G{E-E)*A%e'^^-'^'/^dE, x > a, .^. 

y^j 

i)lAx,t) = J G{E-E)*a'^e^'^-'^'/^dE, < x < a. 



Here, for forming the expression for the WP refiected from the boundary, we select 
the increasing part of stationary solution WF a°exp(^x) only. Imposing the condition of 



'? (x t) = 


= ^Ux,t), 




< X < a, 


iix,t) -- 


+00 

= IG{E- 



- E) * A\e-'^^-''^^/^dE, 


X <0, 


'refi^,t) -- 


+00 
= IGiE- 



- E) * p^e-^'^-'^'/^dE, 


< X < a. 



continuity on the time-dependent WF and its derivative at point x = a, we obtain 2 new 
equations, from which we find the unknowns coefficients Aj^ and a". 

At the third step the WP, tunnehng in region II, is incident upon the first boundary 
of the barrier. Then it transforms into the WP, transmitted through this boundary and 
propagated further in region I, and into the WP, refiected from boundary and tunneled 
back in region II. For WPs one can write: 



(6) 



Using the conditions of continuity for the time- dependent WF and its derivative at 
point X = 0, we obtain the unknowns coefficients A\ and /3^. 

Analyzing further possible processes of transition (and reflection) of WPs through the 
boundaries of barrier, we come to a deduction, that any of the following steps can be 
reduced to one of 3 considered above. For unknown coefficients a*, /3*,74^ and A^, used 
in expressions for WPs, forming in result of some internal refiection from the boundaries, 
one can obtain the recurrence relations. 

Considering the propagation of WP by such way, we obtain the expressions for WF 
on each region which can be written through the series of multiple WPs. We determine 
the resultant expressions for the incident, transmitted and refiected WPs in relation to 
the barrier: 

+00 _ 

iJ^nc{x, t) = J G{E -E)* e'^^'-'^'/^dE, x<0, 



+00 _ +00 

^Pt^{x,t) = J G{E -E)*Y: Afre'''''-'^'/^dE, x > a, (7) 

j=o 

+ 00 _ +CO 

ipref{x,t) = J G{E -E)*Y: A^j^e-'^''-'^'/^dE, x<0. 

i=0 

The series of coefficients a", /?", A!^ and A"^ can be calculated on the basis of obtained 
reccurence relations and have the following form: 

EA^ = -^ , Ea" = ^\ ^ , 

n=0 „ r sub ™=0 rsub 



kl * D_ +^ 2k{ii + k) 



where 



+00 ^ 

n=0 rsub ri=0 P sub 

Fsub = {k^-e)D- + 2tk^D+, 

Kq = k + 1, = 2 ■ 



(9) 



All obtained expressions for the series E«"; J2P"', E^t ^^^ ^A]^ coincide with 
the corresponding coefficients a, (3, At and Ar of the expressions (1), calculated by the 
stationary methods ISEj- Using following substitution 

Ze ^ fc2, (10) 



where k2 = ^j2m{E — Vi) is the wave number for the case of above-barrier energies, 
expression for coefficients a", /9", A^ and A^ for each step, expressions for WF for each 
step, the total WP transforms into the corresponding expressions for the problem of the 
particle propagation above this barrier. Besides the following property is fulfilled: 



+00 



n=0 



2 



+ 



n=0 



1. fill 



3 Tunneling and reflecting times 

Let's find the time of leaving of WP from the barrier formed in result of n refiections from 
the boundaries (we name such WP as n-fold packet). Using the expressions (4) for the 
first step, one can determine the equation describing the propagation of the maximum 
(peak) for incident, transmitted and refiected WPs in relation to the first boundary: 

d d d 

■^arg^pinc{x,t) = -^arg^tr{x,t) = ■^arg^pref{x,t) = const. (12) 

Let the WP in region I be incident upon the first boundary of the barrier in the initial 
time moment tmc- One can find the time tl^f of the leaving of the refiected WP from the 
first boundary to region I: 



-, ^ darqA% 2m 



^ref ~ '^inc + '^ 777=, — tine + TTT' K^'-') 



Similarly, for the time t\^ of the leaving of the transmitted WP from this boundary to 
region II we receive: 

-^ Bar g 13^ m , , 

ttr = tinc + n -Q^ =tinc+J-j^- (14) 

Further we consider the second step, when the WP transits through the second bound- 
ary of the barrier at point x = a. Using the equations (5) and (12), one can obtain the 
time of the leaving of the transmitted WP to region III and refiected WP to region II 
from the second boundary: 



„ ma ^ darqA% 2m 
2 = 7-. -I u h - L = f. _| 

^"^^ hk^ dE ^"""^ hki' 



f^ = f. -\ u fi - L = f. _|_ 

n,k oE 

darga^ 3m, 






Using this approach further, we find the time of the leaving of the n-fold WP from 
the second boundary to region III for the even step n: 

ma ^ daroA^ , ., 2m, 

Similarly, for the odd step n (except for first step) we find the time of an leaving of 
the n-fold WP from the first boundary to region I: 

^daroA^ . ,4m ,_„, 

tref=Unc + n—^^=t,nc+{n+l) — . (17) 

5 



The complete WF, describing a particle transmitted through barrier, represents total 
WP, the coordinates {x, t) of maximum of which can be identified with coordinates of 
a particle (its most probable position), transmitted through the barrier and propagated 
further in region III. The time corresponding to the leaving of this maximum one can 
identify with the phase tunneling time of particle through the barrier. One can obtain: 



Ph _ + + _ "'" I fc n=0 



ma 



hoo 

darg E -4^ 



T 



T-i^n — + _ + — _|_ fe 



Tik qE 



m k^sh{2^a) + 2a^k^{^^ - k^) 



hk^ 4fc2^2 + k^sh^^a) 

The phase time of reflection of a particle from the barrier takes into account only all 
WPs leaving from the barrier in region I: 

+CX1 

darg E ^r 

Ph J. J. 1- n=0 Ph n QA 

Ve/ ~ 're/ '^inc — '' 0771 ~ 'tun- \^^) 

Considering high enough (and wide enough) barrier (limit: ^a ^ 1), one can obtain 
the following expression for the phase times of tunneling and reflection (2]: 

^Ph - ^Ph ]:_ /rjQN 

'tun — 'ref ^ p^ K^^ J 

where v = hk/m is the group velocity. One can see, that these times at such limit 
coincide with time of reflection tl^r at the first step and with time of transition t^^ at the 
second step. Therefore, considering n-fold transmitted and reflected WPs at such limit 
the maximum of total WP is determined by first packets only. 

The phase times of tunneling and reflection, calculated using of the method of multiple 
internal reflections, convergent with results of 0. 

4 The particle tunnels through the ID- double-humb 
rectangular barrier 

The method of multiple internal reflections allows simply enough to calculate the ampli- 
tudes of transmitted and reflected WPs in relation to the barrier for normalized incident 
WP, if the general stationary solutions of WF are known. This method is effective in 
finding the amplitudes for stationary WF. The solutions obtained by this method for 
the problem with barriers of more composite form, than rectangular, appear easier, than 
solutions obtained by stationary methods. 

Let's consider the particle propagating above the one-dimensional double-humb rect- 
angular barrier (see Fig. E}. Solving the stationary Schrodinger equations in each region, 
one can obtain the general solution for stationary WF: 



^{x) 



gjfcx _^ ^^g ikx^ ^Qj^ X < Ri, 

^^^-iksx ^ p^^ik^x^ fo^ R2<x < Rs, (21) 

Q,^Q-ik4X j^ IJ^eik4X^ fQj, i?3 < X < R4, 

Atc''^^, for X > i?4. 



where k = ^y/2mE, ki = ^J2m{E — Vi) is the wave numbers, i is the number of region. 
For finding the unknown coefficients At, Ar, Ui and f3i one can impose the continuity 
condition for WF and its derivative at points x = Ri. In resuh one can obtain 8 equations, 
from which the unknowns coefficients can be calculated. In this the standard stationary 
solution method consists. 

Let's find these coefficients, using the method of multiple internal refiections for so- 
lution of this problem. Analysing 7 undependent steps of WP propagation through the 
barrier boundaries and finding the reccurent relations between unknown koefficients A^, 
A^, a'j and f3j, we can write the expressions for incident, transmitted and refiected WPs in 
relation to the barrier. The sums of coefficients can be calculated on the basis of obtained 
reccurent relations and have the following form: 



[32 

At 



n=0 



n=0 

+ 00 

n=0 

+ 00 

E(32 

n=0 



Ti^ 



I- R^Fri' 



l-R^Fri' 



Fit 



1 - Fr-oFri ' 



FU 



1 - Fr^Ri ' 



n=0 



(22) 



where 



Fr:l 



Fr; 



Fti 



Rn 



Ro + 



1 - FrtR2 ' ^ 
T2R1T+ 



nr^ 



1 — R2 Ri 



+ D-' 



1 - R^Ri 



Fr:, 



Ftt 



Rt+ ^ ^ -^ 



RK + 



1 — /?4 -R3 

T^FroT^ 



2 -^3 



1 - RtFr-n 



FttT^ 



1 - RtFr. 



The coefficients T^ and Rf have the following form: 



(23) 



T+ 



%' 



T' 



n 



2k 



Rt 
Rt 

Ra 



k + k2 



J(k-k2)Ri 



k2 + h 

2k 



ks + kA 
2kA 



i(k—k4,)Rn 



k4, + k 
_ k — k2 ^2ikRi 



k + ko 
k-}. - k 



h + h 



^g2ifc2R2 



ks — k. 

ki — k 
k4 + k 



lg2ifc3R3 

ks + ki 

.j2ik4R4 



T{ 



2k 



J(k2-k3)R2 rp- 



2k2_ 

k 
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T4- 

RI 
? R2 
) R3 

Ra 



2 ^-i(k-k2)Ri 



k + k2 

^^3 ^-i(k2-k3)R2 

k2 + h 

2k4 



ks + ki 

2k 

1 re 

ki + k 

k2-k 



k2 + k 

^3 - ^2 ^-2ik3R2 

kz + k2 

kj — ks 2ik4R3 



-i{k3-k4)Rz 



-i{k—k4)R4 



-2ik2Ri 



ki + /Cs 

k — ki^ 

k + kd 



-2ikRA 



(24) 



(25) 



From these expressions using the substitution similar (10), one can obtain the corre- 
sponding solutions for the problem when the particle tunnels under the barrier having the 
same form. 

Applying the method of multiple internal reflections by this way, one can find the 
expressions for WPs (and also amplitudes for corresponding stationary WF) in the above- 
barrier regions for an one- dimensional problem with a barrier consisting of an arbitrary 
finite number of potential steps. To receive such solutions using the usual stationary 
approach it appears much more complicatedly. This perspective of the method of multiple 
internal reflections is found out for the first time and is displayed the more brightly, if the 
barrier have more composite form. 

5 The scattering of a particle on spherically symmet- 
ric rectangular barrier 

Let's consider the problem of scattering of a particle on the spherically symmetric field, 
the radial part of which looks like (see Fig. EJ: 

( -Vo, for r < Ri 
V{r) = \ Fi, for i?i < r < i?2, (region 11); (26) 



K), 


for r < Ri, 


(region I); 


^1, 


for Ri <r < i?2, 


(region 11); 


0, 


for r > i?2, 


(region III). 



We study the case, when the moment / = and the energy levels lay below than 
the height of barrier. Let's consider the standard stationary method of problem solution. 
Solving the stationary Schrodinger equation in each region, we obtain the general solution 
ofWF: 

^(r,^,V^) = ^*l^^(^,<^), (27) 

r 

^(•g-*fcir _ gjfcir^^ for r < i?i, (region I), 

X{r) = { ae«^ + fSe'^'', for Ri < r < i?2, (region II), (28) 

3-jfcr _|_ g^tkr^ -j?Qj^ r > i?2; (rcgiou III), 



where Yim{0,ip) is the spherical function, ki = jj2m{E + VQ), ^ = |w2m(Vi— £J), 

k = ^y2mE. One can find the unknown coefficients S, A, a and /5 from the continuity 
condition for WF and its derivative at points r = Ri and r = i?2- With the stationary 
point of view the tunneling of particle can be described on the basis of spherical vawes. 
In the obtained solution it is impossible to separate the vawe transmitted through the 
barrier from wave reflected from the barrier, and, therefore, to separate their amplitudes 
(only one item S'e*'^^ contains both transmitted, and reflected waves). 

The method of multiple internal reflections allows to find the solution of this problem. 
Let's apply it to this problem. We study the propagation of WPs sequentially on steps of 
its transition in relation to each of boundaries of the barrier (similarly to one-dimensional 
problem) . 

In result of multiple internal reflections (and transitions) in relation to the boundaries 
of the barrier the total time-dependent WF in each region can be written in form of series. 



composed from the incoming and outcoming WPs. One can calculate the series: 



+00 

n=l 

+00 

E A 



^S h ^2ikiRi\^2({R2-Ri)-2ikR2 



\?4 + h 



Fsub * {k + iO^ 

/^jJ^f^-ikR2+tkiRi-^{R2-Ri) 



n=0 Fsub *ik + iQ {ki + iC) 

ki + i^ 



(29) 



+00 

E «" -- 

n=0 

n=0 Fsub 



Fsub *{k + i^) 



) 



where 



ki + i^ {k + i^){ki + i^) /ggx 

_ ^ ~ ^^ ^-2C(R2-Ri)+2ikiRi 

k + ti 
Now we determine the incident, transmitted and reflected WPs in relation to the 
barrier at whole. Considering them for region III, we write: 



rOO 



Xi 



mc 



r, t) = J g{E - E)e{Vi - E) * e-''''-'^*/^dE, 







Xtr{r, t) = J g{E - E)9{Vi - E) * Str * e'^'-'^'/^dE, (31) 





roo 



where 



X.e/(r, t) = j g{E - E)e{V, - E) * Srefe^'^^'^^'/^'dE, 




+ C)0 

Str = E S"", 



Jo (32) 



Jref 

The expression S represents the S-matrix. Thus, using the method of multiple in- 
ternal reflections it appears possible to divide it into two components corresponding to 
amplitudes of transmitted and reflected WPs in relation to the barrier at whole. This 
property having physical sense, is obtained for the first time. 

Expressions for coefficients S*", A", a" and /5" for each step, expression for WF for each 
step, the series of coefficients 5*", A", a" and /5" under the substitution (10) transform 
into the corresponding expressions for the solution of problem of WP propagation above 
the barrier. The series of coefficients S*", A", a" and /3" coincide with the corresponding 
coefficients S, A, a and P, calculated by stationary methods. 

Similarly to the equations (12) for the one-dimensional problem one can determine the 
equation for propagation of the maximum of incident, transmitted and reflected WPs in 
relation to the barrier for the spherically symmetric problem. One can obtain the times 
necessary for transitting WP through the barrier and for reflecting WP from the barrier: 



Ph _ , , _2mR2 dargStr_ 

hk ^ BE ' 
Ph _ + + _ 2mi?2 , oargSref 



Ttun ^tr l^nc f^^ ^ BE 



(33) 
^Ph. _ 4. _ J- — ZlLHIl I t, ""' 

ref 're/ '^inc j- , ~r 't 



hk BE 



For the problem of particle tunneling under the barrier we receive: 
P,, ^ _d_ i^-h- {i^ + fcQe^^^^^i 






2m 



and for the problem of WP propagating above the barrier we receive: 

p. 2m(i?2 -Ri) ^ d k2-h- {k2 + ki)e^'^^^^ 
nun = TT + ^TTT;ar(7- 



^Ph _ Q 



nk2 dE ^ {k + k2){ki + k2)Fabove ' (35) 



where Fabove can be obtained from the expression for Fgub using the substitution (10). 

Let's consider the enough high and wide barrier (limits: ^{R2—Ri) — > +00, ^ -^ +00). 
For tunneling time we obtain the following expression: 



_p^ 2m , 4mi?isin2A;ii?i(l -2cos2A;ii?i) 
'*™ hk^ ' ;ie(l-cos2A;ii?i; 



<t = ^ + ^"^^" .r""^ :T ""^ (36) 



The tunneling time does not depend on barrier width (the effect of Hartman-Fletcher), 
but depends on ki and Ri. 

6 Conclusions 

At the first time the non- stationary method, describing a process of tunneling of a par- 
ticle through the barrier on the basis of consideration of multiple internal reflections of 
wave packets relatively with the barrier boundaries and named as the method of multiple 
internal reflections, is presented in this paper. In accordence with this method one can 
describe the tunneling process in dependence on time and study the specific features of 
tunneling at any interesting moment of time and point of space in details. 

The one- dimensional stationary problem of tunneling of the particle through the rect- 
angular barrier with taking into account of the multiple internal reflections was earlier 
solved IHE]- The new specific perspective of the method of multiple internal reflections 
in solving of this problem is a possibility to fulfill the time analysis of tunneling and 
reflection in relation to the barrier (with consideration of multiple internal reflections of 
wave packets). 

Using the method of multiple internal reflections the problem of tunneling of a particle 
through spherically symmetric barrier is considered at the first time and splved. Here, 
using this method it is possible (as against the known stationary approaches) to separate 
the wave packet transmitted through the barrier, from the wave packet reflected from 
the barrier (both packets are spherically divergent). In result, one can calculate such 
stationary parameters as the coefficient of penetrability of particle through the barrier 
and the coefficient of reflectivity of particle from the barrier. 

In result for the spherically symmetric problem for S-matrix the following property is 
fulfilled: 

i. e. it consists of two components corresponding to the transmission and reflection of 
particle in relation to barrier. This property has physical sense and is justifled. 
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Fig. 1: An one-dimension one-humb rectangular barrier 
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Fig. 2: An one- dimension double-humb rectangular barrier 
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Fig. 3: A spherically symmetric rectangular barrier 
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